Probing interfacial pair breaking in tunnel junctions based on 
the first and the second harmonics of the Josephson current 
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It will be shown that a pronounced interfacial pair breaking can be identified in Josephson tunnel 
junctions provided the first ji c and the second j2c harmonics of the supercurrent, as well as the 
depairing current in the bulk jd P , axe known. Namely, within the Ginzburg-Landau theory a strong 
interfacial pair breaking results in the relation j2cjd P S> ji c , while in standard junctions, with 
negligibly small pair breaking, the relation of opposite character takes place. 

PACS numbers: 74.50. +r, 74.20. De 



I. INTRODUCTION 

A remarkable property of superconducting weak links 
is that the local conditions in a small transition region 
control the whole process of charge transport. For the 
same reason, interface-induced suppression of the super- 
conducting condensate density can have a considerable 
influence on the Josephson effect. A strong surface pair 
breaking has been theoretically established for various 
unconventional superconductors as well as for magnetic 
interlayers under certain conditions 1-9 . Therefore, prob- 
ing the condensate density near the interface would pro- 
vide valuable information for studying and controlling 
fundamental characteristics of the superconducting junc- 
tions. It is still an ongoing problem for the junctions 
though the order parameter profiles near superconductor- 
vacuum surfaces have been recently determined using a 
Scanning Tunneling Microscopy method with a supercon- 
ducting tip 10 . 

In superconducting tunnel junctions the first harmonic 
ji = jic sin X usually strongly dominates the Josephson 
current j = j\c sin X + J2c sin %X + • • • i while the second 
harmonic j 2 = ]2c sin 2x represents a small correction to 
the first one, \j 2 c\ ^ | ji c | , mostly due to a small junction 
transparency. Qualitatively different phase dependencies 
of the two harmonics allow to study and distinguish be- 
tween them experimentally. It is therefore of interest to 
find out, which characteristic properties of the supercon- 
ducting junctions can be identified with the data pro- 
vided by the two harmonics. Thus the first harmonic, 
as opposed to the second one, is known to be notice- 
ably suppressed both at 0-7T transitions as well as in the 
junctions involving unconventional superconductors with 
special interface-to-crystal orientations 11-18 . However, 
except for these special cases, the relation \j 2c \ <?C |ji c | al- 
ways takes place and does not qualitatively discriminate 
between various superconducting tunnel junctions. 

This paper suggests a test, which will be derived within 
the Ginzburg-Landau (GL) theory and will allow identifi- 
cation of a pronounced interfacial pair breaking in tunnel 
junctions, provided the first and the second harmonics, as 
well as the depairing current in the bulk jd p , are known. 
The relation j 2c jd P ^> j\ c will be shown to take place in 
tunnel junctions with a strong interfacial pair breaking, 



while in standard tunnel junctions, with negligibly small 
pair breaking, it will be j 2 cjdp = 0.27ji c < j\ c . The spe- 
cific temperature dependencies of the two harmonics near 
T c will also be determined. The self-consistency is shown 
to alter existing estimates of both harmonics consider- 
ably. Initially, the theory is based on the interface free 
energy, containing only the terms, which are quadratic or 
bilinear in the superconducting order parameters. Later 
on, the study will be extended to include the next order, 
quartic and biquadratic, corrections. They will be shown 
to result in the material dependent coefficients, which are 
independent of temperature near T c and should, gener- 
ally, be kept on a par with numerical terms of the order 
of unity in the expressions for the order parameters. 



II. BASIC EQUATIONS 

Consider tunnel junctions with the spatially constant 
width, which is much less than the Josephson penetration 
length, and with a plane interlaycr at x — of zero length 
within the GL approach. Assume the usual form of the 
GL free energy, which applies, for example, to s-wave 
and d x 2_ y 2-wsve junctions: T = Th\ + J~b 2 + Fint- Here 
•^61(2) are the bulk free energies of two superconducting 
leads and J-int is the interface free energy. For a junction 
with two identical superconductors, the bulk free energies 
have identical coefficients 

•^61(2) =/ (^|V* 1( 2)| 2 + a|*i( 2 )| 2 + 

J ^1(2) 

+ (6/2) |*i(2) f) dV 1{2) . (1) 

Here K, b, a > 0, a = ar = -a(T c - T)/T c . 

Asymmetry can be generally maintained by different 
conditions on the opposite sides of the interface, as in d- 
wave junctions with different crystal-to-interface orienta- 
tions, and/or in junctions with asymmetric magnetic in- 
terfaces. Then the interface free energy incorporates dif- 
ferent contributions from the two superconducting banks: 
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In addition to the main terms, which are quadratic or 
bilinear in the order-parameter moduli, the quartic and 
biquadratic terms of the next order of smallness near T c 
are kept in (2). In tunnel junctions with small trans- 
parencies D<1 one gets 312, T71, 172 oc V, h%2, fu oc V 2 . 

For the order parameter f(x)e lx ^ normalized to / = 1 
in the bulk without superflow, the first integral of the GL 
equation in the presence of the supercurrent 19 takes the 
form 



l /4+ J^i 
2 J 27 f 2 



(3) 



Here x = £ = \jK/\a\ is the temperature de- 

pendent superconducting coherence length, j is the spa- 
tially constant normalized current density j = j/jdp = 
- (3^3/2) (dx/di)/ 2 , j dp = 8|e|iT 1 /a| |3/2/3V'jffi 6 is the 
depairing current in the bulk, and foe is the asymptotic 
value of / in the depth of the superconducting leads. 

The boundary conditions (BC) originate from the vari- 
ation of Tint and from the bulk gradient terms integrated 
by parts. One starts with the BC in standard linear ap- 
proximation in /i(-0) = /10 or / 2 (+0) = / 20 - Taking 
real and imaginary parts of the BC for the complex quan- 
tity f(x)e lx ( x \ one finds the following linear BC and the 
expression for the supercurrent 

(dfi/dx) = (-If [(gu + <?i 2 )/io - .912 cosx fj ] , (4) 
j = (3^3/2)312/10/20 sin x- (5) 

Here i = 1,2, i = 3 — i. The phase difference of the or- 
der parameters across the interface is x = X10 ~ X20, 
and gi2 = gu£(T)/K, 311 = g n ^(T)/K and 322 = 
g22^(T)/K are the effective dimensionless coefficients. 

For tunnel junctions I312I *C 1- Since the order param- 
eters near pair-breaking interfaces vary on a scale > £(T), 
one gets from (4) guf ia < 1, on account of \g12\fm < 1- 
This signifies, in particular, that for j,i > 1 a strong 
interfacial pair breaking fio < g^ 1 1 occurs. 



III. TEST FOR A PRONOUNCED 
INTERFACIAL PAIR BREAKING 

Consider the supercurrent within the second order per- 
turbation theory in 312. Then, according to (5), quanti- 
ties /10 and /20 should contain the terms of the zeroth 
and the first orders of smallness. One takes x = ±0 in 
(3) and substitutes there (4) and (5). Since the depairing 
effects in the bulk would contribute to (3) only beginning 
with the second order terms, within the given accuracy 
/oo = 1. Then one obtains the following equations for 
f w and / 20 

fa - 2(1 + g\ )fl + 1 - l3,j//„/,., (f i0 - f l0 cos X ) ■ (6) 



In the zeroth order in 312 the solutions arc 

/W 2 = (l/2)( v / 2 T^-3 4l ) 2 , i = l,2. 



(7) 



Eq. (7) involves two solutions of (6). At gu > it de- 
scribes a pair breaking fo < 1, and then |d/(°)/<£c| < 
1/V2. At gu < an enhanced superconductivity at 
the boundary /o > 1 occurs 20-24 . Then the quantity 
\df^'/dx\ Q > V2gfi can take large values, and a strong 
enhancement would induce a characteristic scale substan- 
tially less than £(T) of the leads (see Appendix for de- 
tails) . 

The first order corrections ,/oi ~ /oi + foi \ which 
follow from (6) and (7), are 



(8) 



Substituting the order parameters /oi in (5), one finds 
the first and the second harmonics of the supercurrent 
j = jci sin x+jc2 sin 2% in the Josephson tunnel junctions 



3d- 



3^315 



2 + 5u - .9ii 



1 - 



.912 



2 + 322 ~ 322 

312 



\/2 + Wi V2 + .922 



(9) 
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(10) 



The second harmonic (10) is induced by the proximity 
across the interface. At gu < 0, the quantity |<7jj| is here 
assumed not to be too large to retain \j c \ <C 1 and \gu\ <C 
y/Ka. Otherwise, Eqs. (9) and (10) are applicable at any 
values of gu . Further, the small second and third terms 
in the square brackets in (9) will be neglected. 

One finds from (9) and (10) the following relationship 
between the amplitudes j C 2 and j c \ : 



Jc2 



ill 



E 



9l 



(11) 



Under the conditions \gu\ -C 1, (i = 1,2) one can 
disregard the interfacial proximity effects. Then in the 
original units j cl oc |r|, j c2 oc y/\r\, while the relative 
magnitudes of the two harmonics are described by the 
equalities 



jc2jd P = 0.27j 



cl) 



3c2 = 0.73i2jcl. 



(12) 



Consider now asymmetric junctions with a pronounced 
interfacial pair breaking on one side of the interface, when 
|.9ii I *C 1 and g\ 2 ^> 1, 322 > 0. In d-wave junctions 
this can take place for interface-to-crystal orientations, 
which are close to (100) and (110) orientations on the 
opposite banks of a smooth interface. Then Eqs. (9) 
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and (10) are reduced to j c i « 3v / 3.9i2/(2v / 2.g22), jc2 ~ 
3V3.g? 2 /(45 22 ) while jci oc |t| 3 / 2 , j C 2 oc |r| in the original 
units. The relationships between the harmonics are 

jc2jd P = 0.385<722Jci> 3c2 = 0.7gi 2 jd- (13) 

In symmetric junctions with gn = g 22 > and gf { 3> 1 
one gets from (9) and (10) j cl = 3\/3gi 2 / 'ig 2 2 , j c2 = 
3^3512/4522 ■ Hence, j cl oc t 2 , j c2 oc t 2 , and 

jc2 jd P = 0.77 522 jci , 3c2 = (512/ 522) jcl ■ (14) 

In symmetric junctions the quantity jci/ jcl oc g 22 di- 
minishes with increasing pair breaking 26 . 

It also follows from (9)-(ll) that j c2 jdp = 0.136 j 2 i for 
|.9ii I "C 1 and g 22 S> 1, 522 < 0. In symmetric junctions 
with g 2 t > 1, g u < one gets j c2 j dp = 0.19 j 2 cl j \g 22 f ■ 
Under the conditions 511 < 0, 522 > 0, g u S> 1 (i = 1, 2) 
the relation is j c2 jd P = 0.385 522^1- 

Comparing (12), (13) and (14), as well as the results for 
gu < 0, one can conclude that the quantity |j C 2|jdp/j 2 i 
always exceeds unity, when a pronounced interfacial pair 
breaking g 2 ^> 1, gu > takes place on at least one side 
of the interface. At 0Ag 22 ^> 1, the strong inequality 
bc2|jdp ^ jci emerges as a sure sign of the strong inter- 
facial pair breaking. By contrast, i C 2jdp/ici ' s substan- 
tially less than unity for the negligibly weak pair breaking 
or for the enhanced superconductivity, on both sides of 
the interface. 

Though the specific temperature dependencies, deter- 
mined above for both harmonics at different strengths of 
the pair breaking, could be identified near T Cl there arc no 
striking differences between them. At the same time, the 
power-law temperature dependencies of the harmonics 
jd = jd,o\T\ v< actually drop out of (12) - (14), together 
with the dependencies of the effective coupling constants 
gu = ~gaM~ 1/2 = {9a^lK)\r\- 1/2 = 1,2) and of 
the depairing current jd P = jd P ,o| T | 3 ^ 2 - Hence, (12) - 
(14) are applicable to the "low-temperature" amplitudes 
of the GL theory, and then the relation |j C 2,ob'dp,o 3> Jci 
will be valid, if at least one of 5^0 satisfies the condition 
gu,o 3> 1. The latter condition is more restrictive than 
gu ^ 1. Since the quantities gu incorporate contribu- 
tions from a relatively wide angular interval of quasipar- 
ticle momentum directions, they can be quite large in 
anisotropically paired superconductors near T c , but as a 
rule, decrease substantially when the temperature goes 
down 2 . However, this is generally not the case for pair 
breaking effects induced by magnetic boundaries 8 ' 9 . If 
£te,o < 1 while gu > 1, a crossover from |j c2 |jdp > jci 
close to T c to \ jc2 |idp ^ jci wm show up with decreasing 
temperature, as it is described by (11). 

Assume now gu t o 3> 1. The GL "low-temperature" 
values of the quantities usually exceed their actual val- 
ues at T = by about (2 3) times. This concerns, 
in particular, the depairing current: jdp,o/jdp{T = 0) ~ 
2.6 27 ~ 30 . For the standard Joscphson current 31 , one ob- 
tains jd = jci.oM = 27t 3 T c |t|/7C(3) |e|i?7v near T c and 



j cl (T = 0) = ^Ao/2|e|i? N . Hence Jc i. /jci(T = 0) = 
2.66. Despite the value it would have for analyzing the 
experimental results 10 , there still is no microscopic the- 
ory for the effects of strong interfacial pair breaking in 
a wide temperature range. If, qualitatively, no dramatic 
changes of behavior take place and gu,o 3> 1, the re- 
lation I j C 2 1 Jdp 3> jci could remain valid with decreas- 
ing temperature below the GL domain of applicability 
unless anomalous temperature dependencies, if present, 
come into play, e.g., due to Andreev bound states with 
low energies Eb -C Ao. The temperature dependence of 
(|jc2|jdp/jci) in the whole temperature range is of inter- 
est for further theoretical and experimental studies. 



IV. MICROSCOPIC FORMULA FOR 312 

Microscopic expressions for gi 2 and for gu (i = 1,2) 
can be obtained by comparing the Joscphson currents 
of the GL theory with the corresponding microscopic re- 
sults near T c . Consider here standard symmetric SIS 
tunnel junctions with the negligibly small pair breaking 
\gu\ <C 1. Then the GL expression for the first har- 
monic should coincide with the microscopic Ambcgaokar- 
Baratoff formula 31 near T c : j c i = M e \\ a \9i2 / (fob) = 
7r|A| 2 /(4|e|T c i?7v) . Here Rn is the junction resistance 
in the normal state. Since K = H 2 /4m, \a\ — a\r\ and, 
in the absence of the pair breaking, the BCS gap function 
near T c is |A| 2 = 8tt 2 T c (T c - T)/(7C(3)), one obtains 

~gi2 = 2ir z T c mbt:(T)/{7({?,)e 2 haR N ). (15) 

Eq. (15) can be transformed further with the Gor'kov's 
microscopic formulas for b/a 32 and with the junction re- 
sistance expressed via the averaged transparency R^ 1 = 
e 2 k fT> / Air 2 h. Thus for dirty junctions one obtains 
512 = 0.752?£(T)/£, while for pure junctions <?i 2 = 
3tt 2 Pc:(T)/(14C(3)Co) = 1.762^ (T)/£ - Here I 'is the 
mean free path and £0 = fovf /^T c is the zero-temperature 
coherence length. The quantitative microscopic formulas 
obtained here for <? 12 agree with the earlier estimates 26 . 
In particular, in dirty superconductors the ratio £(T)/l 
can easily reach 100 even at low temperatures. Hence, 
for small and moderate transparencies the quantity gi 2 = 
0.752?£(T)/Z can vary from vanishingly small values in 
the tunneling limit considered in this paper to those well 
exceeding 100 near T c , when a substantial anharmonic 
behavior of the Josephson current takes place 26 . 



V. NEXT ORDER TERMS IN THE CURRENT 

The initial expression (5) for the supercurrent can be 
generalized to include the next order terms, which origi- 
nate from the phase dependent biquadratic contributions 
to (2). The resulting formula is obtained after replacing 
512 -> 512 + fji{\a\l h )fw + m(\a\/b)fi in (5) and adding 
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J/12 =ic2,/sin2x, where 

j c2 j = (3V3/2)(H/6)/ 12 A 2 /! 



(16) 



Here f 12 = fi 2 €(T)/K, ft = raffl/K, i = 1,2. Substi- 
tuting the zeroth order quantities (7) in (16), one obtains 

3c2j = (3V3/8)(|o|/6)/i 2 ( v ^ + & - 5ii) 2 x 

X (v /2 + .92 2 2 -922) 2 - (17) 

Both contributions to the second harmonic (10) and 
(17) are of the second order in transparency oc T> 2 , 
but (17) also contains an additional small parameter 
|r| = (T c - T)/T c , since \a\ = a\r\. This allows to 
disregard (17) in studying the regular problem assumed 
above. However, in a number of specific cases the cou- 
pling constant g\2 can vanish for symmetry reason . 
This concerns, in particular, the asymmetric junction be- 
tween identical dj. 2^2- wave superconductors with exact 
(100) and (110) interfacc-to-crystal orientations on oppo- 
site banks of a smooth plane interlayer 11-15,36-41 . An ad- 
ditional element of the point symmetry inherent in such 
a specific system is the reflection in the xz-plane perpen- 
dicular to the interface. Free energy should be invariant 
under the latter transformation, while the d xy -w&ve order 
parameter on one side of the interface changes its sign 
and the d x 2_ y 2-w&ve order parameter on another side 
keeps its value unchanged. Then the expression contain- 
ing 

|*i -* 2 | in (2) is no longer invariant and, therefore, 
the coefficients 512, fji and 772 should vanish in the case in 
question. By contrast, the term containing \ — in 
(2) remains unchanged under the sign reversal of one of 
the order parameters and, hence, the coefficient /12 can 
maintain its regular value. 

In reality, the first harmonic j c \ remains finite and, 
along with j C 2,/, still represents a substantial part of 
the supercurrent, mainly due to intcrfacial imperfections 
such as faceting, roughness, etc 11-15 . Since \j C 2j\ -C jdp, 
the relation | j c i| < \j C 2,f\ always results in the condition 
Iic2,/bdp 3> ici; which consequently loses its importance 
in the special case of strongly suppressed 512. 



VI. NEXT ORDER TERMS IN THE BC 

Let the parameters gi 2 and gu (i = 1,2) be indepen- 
dent of T near T c . Since gi2,gu oc £(T), then close to T c 
one will get \gi2\ 3> 1 and/or \ga\ ^> 1 due to large values 
of £(T). However, the coupling constants \gi2\ and \gu\ 
can themselves be very small and the temperature range 
with large \gi2\ and/or \gu\ be too narrow. While the 
condition | ^12 1 *C 1, resulting in the tunneling behavior, 
is assumed throughout this paper, the range of variations 
of gu , defined by the strength of interfacial proximity ef- 
fects, is quite wide. It contains, for instance, small values 
of \gu\- For this reason numerical coefficients of the order 
of unity, originating from (3), have been kept in (6)-(ll) 



on an equal footing with g 2 ^ However, the additional 
terms of the next order of smallness, which come from 
the BC, can be comparable with the terms referred to 
above and should generally be taken into account. 

To clarify the point, let's represent the BC schemat- 
ically as (dfi/dx)a ~ -Aj,o + (\a\/b)Ais (i=l,2). Here 
Aifl = Ai,o£(T) I K is linear in the order parameters and 
coincides with the right hand side of (4). The correc- 
tion A\ t i = Ai.i£(T)/K appears in the BC both from 
the quartic and biquadratic terms of the interface free 
energy (2) and from the weak temperature dependence 
of the GL coefficients in Aifl. Therefore, in addition, it 
involves the temperature derivatives of the coefficients. 
As (3) contains (df/dx) 2 , let's consider (dfi/dx) 2 , « 

Af + 2(|a|/6)^4j i o^4i,i- Here it is the crossed prod- 
uct, which is the next order correction to the i-th equa- 
tion for the self-consistent order parameters. Since the 
expression 2(\a\/b)A t0 A t i = 2\a\i 2 (T)A l ^A t i/ (bK 2 ) = 
2AmAu I (bK) depends on temperature solely via the or- 
der parameter amplitudes entering i( ), it results in 
temperature independent coefficients in the equations for 
the order parameters. 

The corresponding terms should, in general, be taken 
into account for the quantitative description of the 
Josephson current. However, as the main contribution 
to (3) from (dfi/dx)^ is quadratic and the correction is 
linear in Aio, for sufficiently large |^4io| t ne correction 
is negligibly small as compared to A 2 . For sufficiently 
small |„4jo| the correction can now also be disregarded as 
compared to the coefficients of the order of unity in (3). 

In tunnel junctions, the basic correction of the given 
origin is described by the crossed product Agahu/ (Kb). 
In particular, in the zeroth approximation in the trans- 
parency the order parameters are 



f (o) 2 

JiO 



1 



9ii 



~2^ 

9u) 



Li 



(18) 



where Li = 1 — (Aguha) / (Kb) and gu, ha > 0. The quan- 
g 2 i jK\a\ are implied here and below to involve 

2Tgi ic (dgii/dT) c , 



tities g\ 
g\ in the expanded form g 2 



q 2 



where weak temperature dependence of gu near T c is 
taken into account in linear in r approximation. 

The linear in gi2 first harmonic is obtained by substi- 
tuting (18) in (5). Calculating also the second harmonic, 
one obtains the modified relation between the second and 
the first harmonics: 



Jc2 



2 



a2 
Jcl 



1 



9u 



9-1 



+ V( 1 + 3S) 



(19) 

In disregarding the term (4guhu)/(Kb), Eqs. (18) and 
(19) are reduced to the previous ones, Eqs. (7) and 
(11). Since the two parameters gu = ga^(T)/K and 
4gnhn/(Kb) are independent of each other, the condi- 
tions \gu\ <C 1 do not generally exclude the special case 
Mgu\hu/(Kb) > 1. Then quadratic in hu corrections can 
also be noticeable. However, for sufficiently small \gu\ the 
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opposite conditions (4\ga\hu/(K b)) 1 / 2 <C 1 occur and al- 
low to disregard all the corresponding terms. 

For large gu the term (4guhu)/(Kb) becomes negligi- 
bly small in (19), when the temperature dependent con- 
dition C(^)fffi ^ 4/i,i/6, in keeping with gu^> 1, is valid. 
Also, for one and the same gu, the right hand side in (19) 
is always larger than that in (11), if (guhu)/(Kb) > 0. 
Therefore, the modified formulas do not alter the main 
statement of this paper. 

In conclusion, a test for identification of a pronounced 
intcrfacial pair breaking in Josephson tunnel junctions 
has been proposed and theoretically verified in this paper, 
based on Eqs. (9)-(ll) and (19) obtained within the self- 
consistent theory of the Josephson current. The main 
statement is that the condition Jc2jdp 3* jci indicates to 
a strong interfacial pair breaking at least on one side of 
the interface, if the first and the second harmonics satisfy 
the conventional relation j C 2 <C | j c i \ ■ 
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Appendix: Order parameter profiles near 
impenetrable boundaries 



The expression for the order parameter in supercon- 
ducting half space with the pair producing surface di- 
rectly follows from (A. 2), since for each /(i), which meets 
(A.l), the functon l/f(x) satisfies the same equation 
(A.l). This results in the following solution 



jpp(x) = coth 



V2 



(A.3) 



for which the condition f pp > 1 holds throughout the half 
space x > 0. 

The order parameter /o, taken on the surface and de- 
scribed by (7), can be alternatively determined by min- 
imizing full free energy (1) and (2) with the solutions 
(A. 2) or (A.3). Explicit integration in (1) with (A. 2) or 
(A.3) results in the part of the bulk free energy modified 
by the boundary. Retaining only the quadratic in the 
order parameter term in the surface free energy (2), one 
finds the full free energy per unit square of an impene- 
trable surface 



T 



g|q| 3 / 2 
V2b 



o - 2/o + ^/ 3 



V2.9/ 2 



(A.4) 



for both solutions. The extremum of (A.4) does result in 
(7) irrespective of the sign of g. The surface suppresses 
the superconducting order parameter at g > 0, while at 
g < the superconductivity is enhanced near the surface. 
It follows from (A.2) and (A.3) 



dx 



1 

V2 



scch 



x 



V2 



(A.5) 



The GL theory allows a detailed description of the spa- 
tial profiles of the order parameters near impenetrable 
boundaries. Here the boundaries, which either suppress 
or enhance the superconductivity in their vicinities, are 
considered jointly. 

In the case in question the supercurrent vanishes and 
foo = 1. Then Eq. (3) reduces to 



dm 

dx 



1 



[1 - fix)] 



(A.l) 



The solution of (A.l), which is relevant to the order 
parameter near a pair breaking surface at x — satisfies 
the condition / < 1 throughout the halfspacc x > and 
takes the form (see, e.g., Ref. 42) 



/pb(i) = tanh 



x + x 

~7T~ 



(A.2) 



The parameter i > together with the associated or- 
der parameter value on the surface should be determined 
from the boundary conditions. 



df P p(x) 
dx 



- — = csch 
V2 



X + x 



(A.6) 



As seen from (A.5), the order parameter (A.2), which 
is suppressed near the boundary, satisfies not only the 
condition f p b(x) < 1, but also the relation 



df P b{x) 



dx 



< 



1 

71' 



(A.7) 



For the order parameter (A.3), which is enhanced 
near the boundary, one gets / pp (i) > 1. According 
to (A.3) and (A.6), the smaller the parameter xq, the 
larger both the order parameter / PPj o and its derivative 
\dJ pp/ ^lo taken on the boundary. A large spatial deriva- 
tive \df P p/ dx\ corresponds to a small characteristic scale 
induced in a superconductor in the vicinity of the surface. 

For g < and >• 1 one gets from (7) and (4) 
/o ~ V~2\g\ and |d//dx| sa \/2g 2 . Hence, the effective 
characteristic scale near the surface is xq ~ £(T)/\g\. For 
the use of the GL theory near the surface one assumes 
xq 3> £o- This results in the condition \g\ <C 1/s/r, i.e., 
\g\ < y/Ka. 
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